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We describe electrophoresis of spherical dielectric particles in a uniformly aligned nematic 
medium with a negative dielectric anisotropy. A spherical particle that orients the liquid crystal 
(LC) perpendicularly to its surface moves under the application of the uniform direct current 
(DC) or alternating current (AC) electric field. The electric field causes no distortions of the 
LC director far away from the sphere. Electrophoresis in the nematic LC shows two types of 
nonlinearity in the velocity vs. field dependence. The velocity component parallel to the applied 
electric field grows linearly with the field, but when the field is high enough, it also shows a cubic 
dependence. The most interesting is the second type of nonlinear electrophoresis that causes the 
sphere to move perpendicularly to the applied field. This perpendicular component of velocity is 
proportional to the square of the field. The effect exists only in a LC and disappears when the 
material is melted into an isotropic fluid. The quadratic effect is caused by the dipolar symmetry 
of director distortions around the sphere and is classified as an LC-enabled electrophoresis 
(LCEEP). The nonlinear electrophoretic mobility of particles in LCEEP offers a rich variety of 
control parameters to design 3D trajectories of particles for microfluidic and optofluidic applications. 
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I. INTRODUCTION 

Controlled manipulation and motion of small particles 
has become a topic of great interest over the last decade. 
Drug delivery, macromolecule separation, display of in- 
formation and colloidal assembly are just a few examples 
of potential applications. Control over the particles can 
be achieved in a variety of ways, including the design of 
self-propelled particles. Historically, however, the most 
popular driving agent of particle motion has been the 
electric field. The most known technique is electro- 
phoresis in which the charged particles are moved by an 
applied DC electric field. There is a growing interest in 
finding mechanisms for particle manipulation that uses 
an AC electric field, as with the latter, it is much easier 
to produce steady flows and to avoid undesirable elec- 
trochemical reactions. Clearly, the AC driving mechan- 
ism should somehow break the classic linear Smoluchow- 
ski relationship between the velocity of particles and the 
electric field E [1-7]. One of the interesting develop- 
ments was a discovery that a broken symmetry of the 
particle can result in a nonlinear electrophoresis, with a 
quadratic dependence of the velocity on the field. The 
phenomenon was described for particles moving in an 
isotropic fluid, first by Murtzovkin [8] and then recently 
by Squires and Bazant [9] who introduced the term in- 
duced charged electrophoresis (ICEP). Experimentally, 



ICEP was discovered by Velevs group for Janus spheres 
(half dielectric, half metallic) moving in water [10]. An- 
other interesting aspect of the electrically driven motion 
is whether the velocity of the transported particle can be 
altered from a trajectory that is determined by the dir- 
ection of the field and polarity of the particles charge. A 
particular case is the reversal of the electroosmotic flow 
of isotropic fluid observed when the frequency of the ap- 
plied AC field changes [7, 11, 12]; the phenomenon is not 
fully understood [13, 14]. Electrically controlled motion 
of particles in an anisotropic fluid (better known as a 
liquid crystal, LC) is clearly more complicated than in 
the isotropic fluid. A replacement of an isotropic fluid 
with a LC should bring first of all some anisotropy in the 
velocity of particles transport parallel and perpendicu- 
lar to the director n (the average direction of molecular 
orientation), associated with the different Stokes drag in 
these two directions, see the review by Stark et al. [15]. 
There are, however, some qualitative differences. Jakli 
et al [16, 17] described both rotational and translational 
motion of colloidal particles in LCs, attributing the ro- 
tational aspects to Quincke mechanism and showing a 
conversion of this rotation into a translational motion as 
a result of hydrodynamic interaction with the walls. Di- 
erking et al [18] has reported on electromigration of mi- 
crospheres in nematic LCs and noticed that the particles 
move under the AC field in the direction perpendicular to 
E and parallel to n; the velocity was measured to have 
a linear dependency on the field. A classic quintessen- 
tial LC effect is the so-called backflow, i.e. generation of 
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LC flow by director reorientations triggered by an elec- 
tric field applied to a LC with dielectric anisotropy. The 
backflow effect was shown to be effective for transport 
and assembly [19, 20] of colloids in a LC; the prime driv- 
ing force of transport is the field-LC interaction rather 
than the field-particle interaction. Ryzhkova, Podgornov 
and Haase [3] have reported on the classic electrophoretic 
experiment staged in a LC in which the electric field was 
acting on charged particles and found a nonlinear (cubic) 
term in the dependence of v on E, in addition to the clas- 
sic linear term. Finally, our group reported that the elec- 
trophoretic motion of a particle in an aligned LC has a 
pronounced nonlinearity with a quadratic dependence of 
v on E [21]. The electrophoretic motion was observed for 
particles that are spherical (an effect that is impossible 
in an isotropic fluid) and for particles that have no elec- 
tric charge. This liquid crystal-enabled electrophoresis 
(LCEEP) was attributed to the dipolar character of dir- 
ector distortions around the particle. The asymmetry of 
the host in LCEEP plays a role similar to the asymmetry 
of particles themselves in ICEP described by Vclcv ct al. 
for Janus spheres moving in water [10]. In our previous 
work on LCEEP [21], we focused mostly on the LC with 
a positive dielectric anisotropy e a = eii — ej_ > 0, where 
ey and ej_ are the dielectric constants measured parallel 
and perpendicular to E, respectively. The electrophor- 
etic motion was observed when the electric field was ap- 
plied parallel to the overall director no (no is fixed in 
space by surface alignment) so that no is not influenced 
far away from the particle, to avoid the backflow-induced 
transport [19, 20]. In this paper, we expand the studies 
to the case e a < 0. This means that the electric field 
does not cause any reorientation as long as E_Lno. Such 
a LC allows us to use two mutually perpendicular field 
directions that do not perturb no far from the particle: 
one perpendicular to the sandwich-type cell, and another 
one parallel to it, Fig.l. We demonstrate that the electric 
field causes two qualitatively different nonlinear electro- 
phoretic effects in the LC. The first effect is a LCEEP 
that drives the spheres along the director no, Fig.l, with 
the velocity growing as vlceep °c E 2 . This effect van- 
ishes when the LC is melted into the isotropic phase by 
increasing the temperature of sample. The second ef- 
fect is reminiscent of a classic nonlinear electrophoresis 
in the isotropic media [22, 23] with the velocity growing 
as i>i3 oc E+aE 3 , where a is a non- vanishing nonlinearity 
coefficient; the nonlinear behavior of velocity is enhanced 
at higher fields and for larger particles. 
Besides the difference in the even vs. odd depend- 
ency of v on E, there is also an important difference 
in the direction of particle motion: the LCEEP velo- 
city is perpendicular to the applied field, vlceep -L E 
(and parallel to no), while the standard nonlinear elec- 
trophoresis velocity is parallel to the field, V13 || E so 
that v LC EEP -L v i3 . 
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Figure 1: Scheme of electrophoretic experiment for spheres 
showing opposite polarities of p y in a LC with e a < aligned 
along rio = (0, 1, 0). (a) Out-of plane configuration with the 
electric field E = (0, 0, E z ); (b) In-plane configuration with 
the electric field E = (^,0,0). The inset shows a polarizing 
microscope texture of a spherical particle with a hyperbolic 
hedgehog near its bottom. 



II. MATERIALS AND METHODS 

We used N, N-didecyl-N-methyl (3-trimethoxysilyl- 
propyl) ammonium chloride (DDMAC) to induce per- 
pendicular alignment of the director n at the particles 
surface. To modify the surface, we dispersed some quant- 
ity of spheres in a solution of DDMAC/Ethanol < 1 wt% 
by sonication during 30 min. The excess of solvent is 
evaporated and the dry-treated spheres are mixed in the 
LC by ultrasonication during 10 min. Polarizing optical 
microscopy textures show that DDMAC yields a normal 
(homeotropic) alignment of LC at the surface of spheres. 
We used dielectric particles of various diameters, namely, 
silica S1O2 spheres of 2a = (5.08 ± 0.5) pm (Bangs Labor- 
atories), borosilicate glass spheres 2a = (9.6 ± l)pm and 
2a = (17.3 ± 1.4)pm (Duke Scientific). They were ad- 
ded to the nematic host in small quantities (less than 1 
wt%) to avoid aggregation. We have chosen the room- 
temperature nematic mixture MLC7026-000 from Merck 
(t n i — 80 °C ) with a negative dielectric anisotropy 
e a = -3.7 (measured at T = 25 °C and 1 kHz). The 
samples represent a thin nematic layer (thickness h var- 
ied from 10 pm to 60p)m sandwiched between two glass 
plates. The glass substrates are coated with polyimide 
PI2555 (Microsystems) which is mechanically buffed to 
produce a uniform alignment along the rubbing direction 
in the plane of the cell, no = (0,1,0) = const, Fig.l. 
The buffing procedure typically results in a small, l°-2° 
pretilt angle between no and the plate. To mitigate the 
possible role of the pretilt in the dynamic phenomena, 
the two plates were assembled with antiparallel rubbing 
directions, which means that no is uniform but slightly 
tilted. The local radial configuration n(r) = f at the im- 
mediate vicinity of the particle needs to match the uni- 
form director no = const away from the sphere. For the 
studied particles, the matching is achieved by a topolo- 



gical point defect, the so-called hyperbolic hedgehog loc- 
ated at short distance from the sphere [24], Fig. 1. There 
is an important symmetry breaking associated with the 
hyperbolic hedgehog, as the director distortions are po- 
lar, Fig. 1. We label the corresponding elastic dipole 
as p and direct it from the hedgehog (negative topolo- 
gical charge -1) towards the sphere (positive topological 
charge 1). The elastic dipole p = (0,p y ,0) is parallel 
to n = (0,1,0) and can adopt two orientations, p y >0 
or p y <0. In a typical sample, the particles with p y >0 
and p y <0 are met with equal probability. Once chosen, 
the sign of elastic dipole cannot change as that would 
require a motion of the hyperbolic hedgehog around the 
sphere, or its splitting into a macroscopic ring, or melt- 
ing the nematic sample into the isotropic phase, with all 
these processes characterized by a huge energy barrier. 
Under a polarizing optical microscope, the silica spheres 
appear as bright circles with a Malthese cross and a 
small tail associated with the hedgehog, Fig. 1 (inset). 
The Malthese cross results from the radial director near 
the sphere. As shown in Fig. 1, the electric field can 
be applied in two geometries with E _L no: (a) out- of 
plane configuration formed by transparent indium tin ox- 
ide (ITO) electrodes deposited onto the glass plates with 
separation between them being 10-60 pm as determined 
by spacers fixing the thickness of the LC slab; (b) in- 
plane configuration formed by two aluminum strip elec- 
trodes separated by 10-15 mm gap. Since e a < 0, there is 
no director reorientations far away from the sphere. The 
in-plane geometry in Fig. 1(b) will be used to explore 
the two different nonlinear mechanisms of electrophor- 
etic propulsion in the LC. The out-of plane geometry in 
Fig. 1(a) is used to demonstrate the ability of LCEEP to 
produce forces that can shift the particles by overcom- 
ing other forces, such as gravity and elastic interaction 
between the director distortions around the particles and 
the bounding walls. 



III. EXPERIMENTAL RESULTS 

(a) Brownian motion. If there is no voltage, the 
spheres in the nematic host experience Brownian motion 
controlled by two different self-diffusion coefficients 
Dn and D±, associated with the motion parallel and 
perpendicular to no, respectively. Loudet et al. de- 
scribed Brownian motion of spheres with quadrupolar 
director distortions [25]. In our case of perpendicular 
surface anchoring and dipolar distortions, presence 
of the topological defect near the sphere breaks the 
" fore — aft" symmetry. However, this feature alone 
does not rectify Brownian diffusion of the sphere and 
does not result in unidirectional movement, Fig. 2. 
The time dependence of the mean square displacements 
(MSD) along and perpendicular to n follows the classic 
linear law, Fig. 2(a), at least for the time scales larger 
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Figure 2: (a) Mean square displacement vs time lag of 2a = 
5.08u.m silica spheres dispersed in E7 LC in the directions 
parallel (y) and perpendicular (x) to the director no measured 
at T = 25 °C; solid lines represent linear fits, (b) Probability 
distribution of particle displacement along and perpendicular 
to the director for time intervals of Is; solid lines represent 
Gaussian fits. 



than 100 ms, with an anisotropic ratio Du/D± = 1.7. 
As these times are larger than the director relaxation 
time, the possible influence of director fluctuations on 
the linear time dependence of MSD vanishes. In the 
absence of E, the time- average displacement of the 
particles averages to zero, as clear from the probability 
distribution P of particle displacements 5 parallel and 
perpendicular to n , Fig. 2(b). 

(b) Dielectric reshaping of the director field near the 
sphere. In Fig. 1, the electric field does not perturb the 
overall director, since e a < 0. Near the sphere, however, 
the director deviates from fio, and is generally not per- 
pendicular to E. The resulting dielectric torque modifies 
the director around the sphere, shrinking the region of 
distortions associated with the hyperbolic hedgehog. 
In principle, this field-induced dielectrically-mediated 
director reorientation can contribute to the propulsion 
of particle. To eliminate this dielectric mechanism of 
propulsion, in our previous work [21], we have performed 
an experiment with a dielectrically compensated mixture 
e a = 0, and found that the particles still move with 



velocity v oc E 2 . The result indicates that besides the 
dielectric reshaping of the director there are other, more 
general, driving mechanisms of LC electrophoresis, for 
example, asymmetric ionic currents around the sphere. 
In the next two sections, we describe electrophoresis in 
a LC for two geometries depicted in Fig.l. 

(c) Electrophoresis in out-of-plane geometry. We 
tracked the position of DDMAC-coated spheres of dia- 
meter 2a — (9.6 ± l)|xm in a well aligned LC, 
fio = (0,1,0), as a function of a vertical DC 
field E z , Fig 1(a), by using fluorescent confocal po- 
larizing microscopy (FCPM) [26] . The nematic LC 
was doped with a very small amount (0.01 wt%) 
of fluorescent dye N,N-Bis(2,5-di-tert-butylphenyl)- 
3,4,9, 10-perylenedicarboximide (BTBP) and the mixture 
was sandwiched between two ITO-coated glass plates 
with a gap h = 50|im between them. The intensity of 
fluorescence is maximum when the transition dipole of 
the dye molecules (parallel to local n) is parallel to the 
polarization of probing light [26]. The silica particles 
appear dark, whereas the hyperbolic hedgehog appears 
as a bright spot next to the sphere, Fig. 3(b), (c). The 
LC environment keeps the colloidal particles in the state 
of levitation [19], since the director distortions around 
the particles are elastically repelled from the bounding 
substrates. The effect facilitates the study of dynamic 
phenomena as it overcomes gravity forces and keeps the 
particles away from the bounding plates. The vertical 
electric field E z displaces the particles along the vertical 
z-axis, Fig. 3, and also along the horizontal y-axis, Fig. 
4. Figure 3(a) shows the field dependence of vertical dis- 
placement (measured between the bottom substrate and 
the center of sphere). There are two mechanisms con- 
tributing to this vertical displacement. First, the elec- 
tric field modifies the director near the spheres and thus 
alters the elastic forces of repulsion from the bounding 
walls [27, 28]. This dielectric effect does not depend 
on the field polarity. The second effect is the polarity- 
dependent electrophoretic shift, Fig. 3(a). TheDDMAC- 
trcatcd spheres move towards the anode, suggesting that 
their electric charge is positive. The horizontal motion 
of particles, caused by the vertical DC field, Fig. 1(a), is 
directed along n and p, i.e. the sphere leads the motion. 
The spheres with p y >0 and p y <0 move in antiparallel dir- 
ections. When the field polarity is reversed, the particles 
continue to move in the same direction, dictated by p. 
This feature indicates that the electrophoretic velocity v y 
does not depend on the surface charge. Furthermore, it 
immediately offers an opportunity to drive the particles 
by an AC field. Figure 4 shows how the velocity v y de- 
pends on the frequency / of the applied electric field for 
cells of different thickness and with different directions of 
p. The dependence v y (f) is non- monotonous with a max- 
imum at / w (0.1 — 0.5)Hz. At low frequencies, v y is par- 
allel to p. As the frequency increases to / « (10 — 50)Hz, 



the particles reverse the direction of motion, so that v y 
and p are antiparallel to each other. However, the motion 
in the reversed direction is slow, Fig. 4, and eventually 
vanishes at / > 100Hz. 
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Figure 3: Electrically controlled levitation of la — 9.6um 
spherical particles in MLC7026-000. (a) z-position of micro- 
spheres in the bulk as a function of the electric field. FCPM 
textures of a vertical cross-section of the nematic cell with a 
colloidal sphere, for (b) negative and (c) positive E z . The 
white scale bar corresponds to 5nm. 

We noticed a small pretilt effect in thick samples: 
particles with p parallel to the rubbing direction at the 
closest (bottom) substrate moved at a slightly higher 
electrophoretic velocity than their counterparts with an 
opposite p. In thick samples, the levitating spheres of 
mass density higher than that of the LC, are located 
closer to the bottom plate [19]. The difference in elastic 
repulsion from the wall experienced by the tilted elastic 




120mV/|im 



*y 



-o 



p <0h=ll.l nm 



n <0h=50.1 urn 
mi I i 



Frequency/(Hz) 

Figure 4: Velocity (y-component) of DDMAC-coated spher- 
ical particles with 2a = 5.08um in MLC7026-000 vs fre- 
quency of the sinusoidal AC field; out-of plane geometry, 
E = (0, 0, Ez) and no = (0, 1, 0). The solid lines are guides to 
the eye. 



dipoles might be the reason for the different electro- 
phoretic mobilities of spheres with the opposite ps. The 
difference vanishes when the gap between substrates is 
reduced (see the data for h = ll|xm in Fig. 4). In thin 
cells, the particles are close to the middle plane (since 
the elastic repulsion from both plates is much stronger 
than the gravity force) [19, 28], and the competing 
effects of the rubbing directions at the opposite plates 
cancel each other. The experiments suggest that there 
are two different mechanisms of electrophoresis in a 
LC, one related to the electric surface charge of the 
particles and causing their motion along the applied 
field, and another one independent of the charge, causing 
the particles to move in the direction perpendicular to 
E. In order to explore these mechanisms in a greater 
detail, we use a different geometry with the in-plane 
electric field E x , Fig. 1(b). The advantage is that the 
electrophoretic forces are not competing with gravity 
and elasticity-mediated forces; the in-plane field E x 
causes electrophoretic motion in the horizontal xy plane 
of the cell. 

(d) Electrophoresis in in-plane geometry. The field 
E = (E X ,Q, 0) was applied using two parallel aluminum 
strip electrodes separated by a gap L = 10 — 15mm. 
For quantitative characterization of electrophoresis, one 
should take into account the finite times of electric char- 
ging in the system. Since LCs always contain some 
amount of ions (caused, e.g., by degradation of mo- 
lecules), impurities from the alignment layers, etc. [29], 
the voltage profile across the cell is time-dependent. 
Once the DC field is applied, the ions move towards 
the electrodes and build electric double layers near them 
(screening the applied field), within a characteristic time 
i~e = XdL/2D [30, 31]. For a typical Debye screening 



length in LC Ad = 0.5u\m [32], L = 10mm, and ionic 
diffusion coefficient D = 3 x 10~ n m 2 /s [33], one estim- 
ates r e = 10 2 s . In the DC field experiment, to avoid 
the electrode screening effects, we measured the particles 
velocities within a fixed time interval At, equal about 3 
min. at low fields (| E |< 10mV/|J.m) and 1 min. or less 
at high fields (| E > lOrnV/um). We experimentally 
established that within the specified values of At, the 
spheres moved with the time- independent velocity. After 
the velocity was measured for a given amplitude and po- 
larity of E, the polarity of the field was reversed and the 
measurement was repeated, within the same time inter- 
val At. The horizontal DC field E = (E x ,0, 0) causes the 
particles to move along the x and y axes, Fig. 5. Con- 
sider first the velocity component v x that is parallel to the 
applied field E x . The dependence v x {E x ) is practically 
linear for small fields, similarly to the electrophoresis in 
regular fluids, described by the Smoluchowskis formula: 



v = /iiE, 



(1) 



where /ii = e m C,/rj is the electrophoretic mobility, C, is the 
zeta potential of the particle, e m and -q are the dielectric 
permittivity and the effective viscosity of the medium, 
respectively. For stronger fields, one clearly observes a 
non-linear behavior with a cubic term: 
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H3E X , 



(2) 



where [13 is a field- independent coefficient. The effect is 

stronger for larger particles, Fig. 5. The cubic nonlin- 

earity for a nematic LC was also reported by Ryzhkova 

et al [3]. 

The velocity component in the direction perpendicular to 

the electric field shows an absolutely different quadratic 

behavior: 



vlceep 



(3 E 2 X , 



(3) 



where (3 is the field-independent coefficient. The sign 
of f3 changes with the reversal of the elastic dipole p. 
The latter feature indicates that the electrophoretic ve- 
locity v y perpendicular to E is not related to the net 
charge of the particles and that the mechanism of this 
quadratic electrophoresis is rooted in the dipolar distor- 
tions of LC around the particle. Once the LC is heated 
into an isotropic melt, the quadratic effect disappears, 
vlceep = v y — 0, while V13 = v x remains nonzero. 

The experiments above suggest that in the nematic 
LC, the electrophoretic velocity- field relationship v(E) 
is tensorial, 

Vi = HiijEj + fajkEjEk + {i 3i jkiEjEkEi, (4) 

with tensorial coefficients /3ijk (rather than a single con- 
stant j3) that relate the velocity components Vi to the 
field components Ej(i,j = x,y,z). Using the experi- 
mental data in Fig. 5 and Eq. (4), we find that for 
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Figure 5: Electrophoretic velocity of positively charged 
DDMAC-coated dielectric spheres of diameter 2a = 5.08(J.m 
and 2a = 17.3um with p y > in the nematic phase (25 °C) of 
MLC7026-000 as a function of the uniform DC electric field; 
in-plane geometry with E = (1^,0,0) and no = (0,1,0); 
(a) z-component of the velocity, parallel to E = (£^,0,0) 
and perpendicular to no = (0,1,0), showing a cubic de- 
pendence on E; (b) y-component of the velocity, parallel to 
no = (0,1,0), showing a quadratic dependence on E. Solid 
lines represent the fitting of experimental data using Eq. (2) 
and Eq. (3) for parts (a) and (b), respectively. 



when a is much larger than the Debye screening length 
Xd [34]. The typical values of Xd in LCs are less than 
1 iim, thus we expect a/ Xd > 1. The second reason 
is LC-enablecl levitation. As shown experimentally by 
Pishnyak et al [35], larger glass particles levitate closer 
to the midplane of the cell than the smaller particles 
(since the elastic forces grow with a faster than the grav- 
ity force). Since the viscous drag on the sphere depends 
on the distance to the bounding plates, the large sphere 
might experience a weaker effective viscous drag than the 
small spheres. Note that in our experiment only some of 
the possible mobility coefficients in Eq.(4) were probed. 
The main reason is that the finite dielectric anisotropy 
of the LC restricts the number of configurations in which 
the field does not perturb the director. By extending the 
studies to the LC with e a = 0, one can explore a much 
wider set of mobility coefficients. A LC with e a — can 
be prepared by mixing two LC materials with e a < 
and e a > in an appropriate proportion, or by using a 
so-called dual- frequency material, in which e a changes its 
sign at a particular frequency of the applied electric field. 



IV. DISCUSSION 

Our experiments demonstrate the existence of two 
types of electrophoretic motion in a LC of negative 
dielectric anisotropy, with two very different velocity- 
field dependences, namely, (1) velocity is an odd function 
of E, viz oc E + aE 3 and (2) an even function of E, 
vlceep oc E 2 . The first effect is known for isotropic 
fluids as Stotz-Wein effect [22, 36]. The second effect is 
observed only if the material is in the LC state; this is 
why we call it an LC-enabled electrophoresis (LCEEP). 
Below we discuss both effects for the geometry of 
Fig. 1(b), in which the experimental velocities follow 
Eq.(2) for the LC analog of Stotz-Wein effect and Eq.(3) 
for the genuine LCEEP, as schematically summarized in 
Fig.6. 



DDMAC-coated positively charged glass spheres of dia- 
meter 2a = 5.08um in the LC MLC7026-000, fi lxx = 
0.011u.m 2 /mV-s, /3 yxx = 1.5 x Kr 4 u.m 3 /mV 2 -s and 
Hzxxxx = 5.2 x 10 _6 |xm 4 /mV 3 -s. For larger spheres 
of the same type, 2a — 17.3u-m, we find [i\ xx = 
0.038^m 2 /mV-s, /3 yxx = 7.4 x 10- 4 um 3 /niV 2 -s an d 
H3xxxx = 1-1 x 10~ 4 |J.m 4 /mV 3 -s. The values of fii 
and /3 are close to those reported previously for pos- 
itively charged borosilicate spheres of diameter 2a — 
17.3|J.m placed in the LC mixture E7 with positive dielec- 
tric anisotropy {^i vv = 0.07|^m 2 /mV-s , (3 yyy = 5 x 
10 _3 |o.m 3 /mV 2 -s) [21]. The data above show that [i\ xx 
tends to increase with the radius a of the spheres. The 
result might be related to two effects. First, within the 
classical theory of electrophoresis in an isotropic fluid, the 
mobility is expected to grow with the size of the particle 



(a) Stotz-Wein effect in a LC. In the experiment, 
the velocity vis = v x , Eq.(2) and Fig. 5(a), follows a 
classical linear dependence on the field, when the field 
is weak. The direction of particle motion changes to 
the opposite when the field polarity is reversed. This 
behavior indicates that the electrophoretic mobility is 
determined by the surface charge of the particle. At 
higher fields, the cubic term in dependency becomes 
apparent, Fig. 5(a). In the isotropic fluids, this Stotz- 
Wien effect [22, 36] is explained by the field-induced 
polarization of the electric double layers that shifts the 
electric potential at the particle surface by a quantity 
oc E 2 , so that V13 — fi±E + fi3E 3 [37]. A similar effect is 
natural for a LC medium, where it can also be enhanced 
by the director realignment in the field caused by the 
dielectric anisotropy, with the realigning torques oc E 2 . 



The field-induced change oc E 2 in the mobility fX\ is 
equivalent to the appearance of the nonlinear coefficient 
^ 3 in Eq.(2). 

(b) Liquid crystal- enabled electrophoresis. The 
LCEEP driven by AC and DC electric fields is observed 
for spherical particles when the director distortions 
around them are of dipolar type, Fig. 6. If the particles 
have a quadrupolar symmetry, LCEEP vanishes [21]. 
These results indicate that the mechanism for LCEEP is 
rooted in the type of director distortions which violate 
the fore-aft (left-right) symmetry of the LC medium 
around the sphere. As already mentioned, anisotropy of 
both dielectric and conductive properties of the nematic 
LC can contribute to the propulsion of the particles 
in the electric field. The dielectric anisotropy of the 
LC causes director realignment in the applied electric 
field that is especially pronounced in the region of the 
hedgehog (as compared to the opposite side of the sphere 
where the director is more homogeneous). Although 
further exploration of this effect is important, we note 
that if the dielectric anisotropy is reduced to zero, the 
LCEEP still persists [21], suggesting that the dielectric 
anisotropy is not a necessary condition for LCEEP. A 
more general mechanism is associated with the different 
mobility of electric ions around the sphere with the 
dipolar director distortions [21]. Consider a dielectric 
sphere with a positive surface charge, screened by a 
cloud of negative counterions, Fig. 6. Besides these two 
types of charges, the LC bulk also contains mobile ions 
of both positive and negative charge (the system as a 
whole is electrically neutral). Once the electric field is 
applied, the mobile ions in the LC start to move in two 
opposite directions, from the two poles of the sphere 
towards its equator. The flow pattern (from the poles 
towards the equator) does not change with the field 
reversal. Since the mobility of ions is anisotropic in a LC 
(mobility along the director is typically higher than in 
the perpendicular direction), the velocity of ions in the 
hyperbolic hedgehog region is expected to be different 
from the velocity at the opposite, hedgehog-free side 
of the sphere, resulting in a broken mirror symmetry 
with respect to the equatorial plane of the sphere that 
is perpendicular to the elastic dipole p. Furthermore, 
since the viscosity of the LC medium is also anisotropic, 
the LC flow triggered around the sphere by ions is also 
expected to follow a broken symmetry pattern, Fig. 6. 

These asymmetries give rise to an electrophoretic ve- 
locity that is polarity- independent, Vlceep = v y oc E 2 . 
For particles moving perpendicular to the electric field, 
the director field asymmetry is the only reason for electro- 
phoretic propulsion since v y is independent of the particle 
charge. Of course, if the particle is charged, then there 
is also a regular electrophoretic effect shown in Fig. 5(a) 
with the velocity V13 = v x parallel (or antiparallel, de- 
pending on the sign of charge) to the direction of the 




Figure 6: Schematic of two types of electrophoretic motion 
of a sphere with normal anchoring at the surface, in a uni- 
formly aligned nematic LC: (a) in-plane (xy) two-dimensional 
representation; (b) three-dimensional view. 



applied field E x , as in Eq. (2). For uncharged particles 
v x = 0, and only the quadratic contribution to LCEEP 
would survive. If the driving field is of an AC type, then 
the average value of Vis = v x vanishes but the component 
vlceep oc E 2 remains non-zero, allowing an AC control 
of LCEEP. The frequency dependence of electrophoretic 
velocity has been extensively studied for ICEP in iso- 
tropic fluids [4, 9, 38] and recently by us for LCEEP [21]. 
In both cases, there are two important time scales: 1) a 
characteristic charging time r c = e m X 2 D / tdD of a dielec- 
tric sphere; 2) the characteristic electrode charging time 
T e = \dL/2D. If only these two scales are relevant, and 
if the Debye layer around the sphere is thin and uniform, 
the velocity shows a frequency dependence 



v(u) 



'(l + u, 2 T c 2 )(l + u; 2 Tf) 



(5) 



with a single peak at the characteristic charging fre- 
quency lo c = t^ 1 of the electric double layer. This 
simplified model cannot capture all the features of the 
frequency dependence of velocity presented in Fig. 4 
which demonstrates two peaks. The model can be 
improved by taking into account that the cloud of 
counterions around the spherical particle in a LC matrix 



is not spherically symmetric; this is a subject of further 
studies. To conclude, we demonstrated a rich variety 
of nonlinear electrophoretic effects when the moving 
particle is located in a LC rather that in an isotropic 
fluid. The electric field is applied in such a way that it 
does not distort the director far away from the particle. 
The electric field causes electrophoresis with velocity 
components that are parallel and orthogonal to the 
field. In the first case, the velocity is proportional to 
the linear and cubic power of the field. In the second 
case, the velocity is quadratic in the field, demonstrating 
a new mechanism of electrophoresis that we call a 
liquid crystal-enabled electrophoresis, or LCEEP. The 
LCEEP occurs when the director distortions around 
the inclusion show dipolar symmetry. The LCEEP 
effect makes it possible to move the particles by an 
AC electric field even when the field is absolutely 
balanced, i.e., not only when its time average is zero 
(E) = 0, but also when the time average of higher 
moments are zero, e.g., (E 3 ) — 0. This represents 
an important difference with the so-called aperiodic 
nonlinear electrophoresis in isotropic fluids with velocity 
that is an odd function of the electric field [39] . It also 
allows one to transport particles of zero charge and 
particles of an absolutely symmetric shape and physical 
properties. The experiments presented above clearly 
demonstrate the tensorial relationship between the 
electrophoretic velocity in the LC and the electric field. 
The fact that the particles move perpendicular to E in 
LCEEP means that the driving voltage can be applied 
across the thickness of the channel (few pm) rather 
than along the pathway (several cm) in electrophoretic 
devices; this effect allows one to use modest voltage 
sources to achieve high fields, which might be especially 
beneficial for portable devices. The possibility of moving 
the particle in different directions without altering the 
direction of the field is also remarkable. The field and 
frequency dependencies shown in Fig. 4 and 5 indicate 
that one can combine two differently oriented driving 
fields with different frequencies to control the overall 3D 
trajectory of particles. Further diversification can be 
achieved by using LCs with pre-pattcrned director field 
and by using LCs of zero dielectric anisotropy e a = so 
that the electric field does not cause the director reori- 
entation regardless of the mutual orientation of fig and E. 
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